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^ ■ Abstract 

We describe a fault-tolerant version of the one-way quantum computer using a cluster state 
in three spatial dimensions. Topologically protected quantum gates are realized by choosing 
appropriate boundary conditions on the cluster. We provide equivalence transformations for 
these boundary conditions that can be used to simplify fault-tolerant circuits and to derive 
circuit identities in a topological manner. The spatial dimensionality of the scheme can be 
! reduced to two by converting one spatial axis of the cluster into time. The error threshold is 

0.75% for each source in an error model with preparation, gate, storage and measurement errors. 
The operational overhead is poly-logarithmic in the circuit size. 



m 

^ ! 1 Introduction 

O . 

The threshold theorem for fault-tolerant quantum computation [U EJ El 0] has established the fact 
Oh that large quantum computations can be performed with arbitrary accuracy, provided that the error 

■y ' level of the elementary components of the quantum computer is below a certain threshold. It now 

becomes important to devise methods for error correction which yield a high threshold, are robust 
against variations of the error model, and can be implemented with small operational overhead. 
An additional desideratum is a simple architecture for the quantum computer, such as requiring no 
long-range interaction. 

The one-way quantum computer provides a method to do this [U [6] , which we describe in detail 
below. We obtain an error threshold estimate of 0.75% for each source in an error model with 
preparation, gate, storage and measurement errors, with a poly-logarithmic multiplicative overhead 
in the circuit size O (~ In 3 $7). It shall be noted that we achieve this threshold in a 2-dimensional 
geometry, only requiring nearest-neighbor translation-invariant Ising interaction. This is relevant 
for experimental realizations based on matter qubits such as cold atoms in optical lattices [7J [8] 
and two-dimensional ion traps [9|, or stationary qubits in quantum dot systems [10] and arrays of 
superconducting qubits [TT]. Geometric constraints are no major concern for fault-tolerant quantum 
computation with photonic qubits [12l E3] . 

The highest known threshold estimate, for a setting without geometric constrains, is 3 x 10~ 2 
|14j . Fault-tolerance is more difficult to achieve in architectures where each qubit can only interact 
with other qubits in its immediate neighborhood. A recent fault-tolerance threshold for a two- 
dimensional lattice of qubits with only local and nearest-neighbor gates is 1.9 x 10~ 5 [15]. We note 
that since the initial work of [16] a number of distinct approaches to topological fault-tolerance 
emerging in lattice systems are being pursued; See [TT] [T8j [T9" l [20] . 

The key element of our method is based on topological tools that become available when the 
dimensionality of the cluster is increased from two to three. In 3D, we combine the universality 
already found in 2D cluster states [21] with the topological error-correcting capability of Kitaev's 
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Figure 1: Topologically protected gates. 



toric code [16]. Then, a one-dimensional sub-structure of the cluster is "carved out" by performing 
local Z-measurements. This leaves us with a non-trivial cluster topology in which a fault-tolerant 
quantum circuit is embedded. Fig. Q] displays topologically protected gates that can be constructed 
in this manner. 

This paper is organized as follows. In the remainder of this section we introduce the necessary 
terminology for discussion of the fault-tolerant QCq- In Section [2] we answer the question of "Why 
can we perform non-abelian gates with surface codes?" and describe topological transformation 
rules on the cluster. We subsequently use these rules to simplify topological circuits. In Section [3] 
we complete the universal set of gates. In Section H] we describe the mapping from a 3D cluster 
state to a two-dimensional physical system plus time. Sections [5] and [6] address the fault-tolerance 
threshold and overhead, respectively. We conclude with a summary and outlook in Section [71 

Before we can start our discussion of the fault-tolerance properties of the QCc, we need to intro- 
duce some necessary notation. This has been done before in [3 [6]. We include a short introduction 
here to make our presentation self-contained. Consider a cluster state \(j))c ° n a lattice C with 
elementary cell as displayed in Fig. [2^,. Qubits are located at the center of faces and edges of C 
The lattice C is subdivided into three regions V, D and S. Each region has its purpose, shape 
and specific measurement basis for its qubits. The qubits in V are measured in the X-basis, the 
qubits in D in the Z-basis, and the qubits in S in either of the eigenbases (X ±Y)/ y/2. V fills up 
most of the cluster. D is composed of thick line-like structures, named defects. S is composed of 
well-separated qubit locations interspersed among the defects. As described in greater detail below, 
the cluster region V provides topological error correction, while regions D and S specify the Clifford 
and non-Clifford parts of a quantum algorithm, respectively. 

We can break up this measurement pattern into gate simulations by establishing the following 
correspondence: quantum gates <-> quantum correlations <-> surfaces. The first part of this cor- 
respondence has been established in [22j . For the second part homology comes into play. The 
correlations of \<j>)c, i.e., the stabilizers, can be identified with 2-chains (surfaces) in C, while errors 
map to 1-chains (lines). Homological equivalence of the chains implies physical equivalence of the 
corresponding operators [S]. This correspondence is key to the presented scheme. Gates are speci- 
fied by a set of surfaces with input and output boundaries, and syndrome measurements correspond 
to closed surfaces (having no boundary). 

C is regarded as a chain complex, C = {C3, C2, Ci, Co}. It has a dual C = {C3, C2, C±, Co} 
whose cubes C3 G C3 map to sites Co G Co of £, whose faces C2 G C2 map to edges c\ G C\ of C, etc. 
The chains have coefficients in Z2. One may switch back and forth between C and £ by a duality 
transformation *( ). C and C are each equipped with a boundary map d, where d o d = 0. 

Operators may be associated with chains as follows. Suppose that for each qubit location a in 
a chain c, a G {c}, there exists an operator E a , and [S a , £&] = for all a, b G {c}. Then, we define 
S(c) := riae{c} Cluster state correlations (i.e. stabilizers) are associated with 2-chains. For 
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the considered lattice, all elements in the cluster state stabilizer take the form K(c 2 )K(c 2 ) with 
c 2 G C 2 , c 2 G C 2 , and 

^(c 2 ) =X(c 2 )Z(0c 2 ), ^(c 2 ) = X(c 2 )Z(cfe 2 ). (1) 

Only those stabilizer elements compatible with the local measurement scheme are useful for infor- 
mation processing. In particular, they need to commute with the measurements in V and D, 

[K(c 2 )K(c 2 ),X a ] = 0, aeV, 

[K(c 2 )K(c 2 ),Z b ] = 0, b£D. { ) 

Due to the presence of a primal lattice L and a dual lattice C, it is convenient to subdivide the 
sets V and D into primal and dual subsets. Specifically, V = V p U Vd, with V p C {C 2 }, Vd C {C 2 }, 
and D = D p U Dd, with Z) p C {Ci}, Dd C {Ci}. With these notions introduced the compatibility 
condition ([2]) may be expressed directly in terms of the chains c 2 and c 2 . If K(c 2 ) and K(c 2 ) have 
support in V U -D only then Eq. ([2|) is equivalent to 

{dc 2 }cD p , {dc 2 }GD d . (3) 



The QCc and surface codes. We need to specify the encoding of logical qubits before explaining 
the encoded gates. For this purpose, let us single out one spatial direction on the cluster as 
'simulating time.' The perpendicular 2D slices provide space for a quantum code. The code which 
fills this plane after the mapping of the three dimensional lattice C onto a 2+1 dimensional one is 
the surface code [23]. 

The number of qubits which can be encoded in such a code depends solely on the surface topology. 
Here we consider a plane with pairs of either electric or magnetic holes. See Fig. [2b. A magnetic 
hole is a plaquette / where the associated stabilizer generator S\j(f) = Z(df) is not enforced on 
the code space, and an electric hole is a site s where the associated stabilizer S+(s) = X(d &s) is 
not enforced on the code space ("#" denotes the duality transformation in 2D). Each hole is the 
intersection of a defect strand with a constant-time slice. 

A pair of holes supports a qubit. For a pair of magnetic holes /, /', the encoded spin flip operator 
is X = X(ci), with {dc\} = {*/,*/'}, and the encoded phase flip operator is Z = Z(c\), with 
ci = df or c\ = df. The operator Z(df + df) is in the code stabilizer 5, 

Z(df + Of) G S. (4) 

For a pair of electric holes s,s' we have X e = X(c'i), with = d^s, Z e = Z(ci), with {dc±} = 
{s, s'}, and 

X (d*s + d*s') G S. (5) 



The simplest gate. Here we illustrate the relation between quantum gates, quantum correlations 
and correlation surfaces (2-chains). We choose the simplest possible example: the identity gate. 

The identity operation is realized by two parallel strands of defect of the same type. We consider 
a block shaped cluster C C C for the support of the identity gate. One of the spatial directions 
on the cluster is singled out as 'simulated time.' The two perpendicular slices of the cluster at 
the 'earliest' and 'latest' times represent the code surfaces I and O for the encoded qubit, with 
1,0 C {Ci} being an integer number of elementary cells apart. As before, we ask which cluster 
state correlations K(c 2 ), K(c 2 ) are compatible with the local measurements in C\(IUO). With the 
additional regions / and O present, the condition ([2]) turns into 

{c 2 } C V p , {dc 2 } C D p UlUO, () 

{c 2 } C VdUlUO, {dc 2 } C D d . W 
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Figure 2: Lattice definitions, a) Elementary cell of the cluster lattice C 1-chains of C (dashed 
lines), and graph edges (solid lines), b) A pair of electric ("e") or magnetic ("m") holes in the code 

6 / 771 G / 7TL 

plane each support an encoded qubit. Z and X denote the encoded Pauli operators Z and 
X, respectively. 




Dp : primal defect 

Figure 3: Correlation surfaces a 2 , a 2 for the identity gate on a primal qubit. 



Surfaces of primal correlations compatible with the local measurements in C\(I U O) can stretch 
through the cluster region V and end in the primal defects and the input- and output regions. They 
cannot end in a dual defect. Surfaces of dual correlations can stretch through V, I and O, and end 
in dual defects. They cannot end in primal defect^. 

We now consider the identity gate on the primal qubit, mediated by a pair of primal defects. The 
relevant primal and dual correlation surfaces are displayed in Fig. O and we denote these special 
surfaces by a 2 and a 2 . Before the local measurement of the qubits in C\(lL)0) the cluster state \<p)c 
obeys K{o 2 )\<j>)c = K{a 2 )\(f>}c = \4>)c- Note that K(a 2 )\ IuQ = Zt^Zo and K(a 2 )\ Iu0 = Xi®X - 
The " 7 " refers to encoding with the surface code displayed in Fig. [2] Thus, for the state \ip)iuo after 
the measurements in C\(I UO), Zj® Zo\ip)iuO — M^iuo and Xj (8) Xo\ip)iuo = =t|V')/uo- This 
is the connection between surfaces (2-chains) and quantum correlations. The connection between 
quantum correlations and gate operation has already been established in Theorem 1 of [22], from 
which the identity gate follows. 

The other Clifford gates (or more precisely, CSS-gates) are derived in a similar manner, invoking 
more complicated correlation surfaces. 

The asymmetry between primal and dual 2-chains in Eq. © arises because I and O axe chosen subsets of C\. 
Physically speaking, we choose the sub-cluster C such that it consists of intact cells of the primal lattice £ at the front 
and back. The cells of the dual lattice L are then cut in half at the front and back of C. 
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2 Topological considerations 



2.1 Why can we perform non-abelian gates with surface codes? 

A limitation of the surface code |16[ [23] is that only an abelian group of gates can be implemented 
fault-tolerantly by braiding operations |16j . In the fault-tolerant QCc, arbitrary CNOT-gates can 
be performed which are non-commuting. Yet, the fault-tolerance of the QCq is based on surface 
codes. How does this fit together? 

The reason why we can do non-abelian gates with a surface code in the QCq is that we change 
the topology of the code surface with time. The preparation of a primal |0) state (dual |+) state) 
introduces a pair of primal (dual) holes into the code surface. The corresponding measurements 
remove pairs of holes. 

The emergence of non-abelian gates through changes in the surface topology can be easily verified 
in the circuit model. Consider first the monodromy of a primal and a dual hole as the means to 
entangle two qubits of opposite type, 



This operation does not change the topology of the code surface. It can be checked with the methods 
described in Section [1] that operation ([7]) acts as a CNOT on the two involved qubits. However, the 
primal qubit is always the target and the dual qubit the control. These gates are still abelian. 

We now supplement these unitary commuting gates with non- unitary operations, namely X- 
and Z-preparations and measurements. They are obviously non-commuting and change the surface 
topology (See Fig. [lb). Can we construct non-commuting unitary operations out of this gate set? 
To this end, we assemble preparations, measurements and the monodromy operation ([7]) to the 
topological circuit displayed in Fig. ®l. It is a deformed version of the gate in Fig. []Jt. Also, it can 
be verified directly in the circuit model that it represents a CNOT-gate (c.f. Fig. [4b). The direction 
of the CNOT can now be chosen freely, and we obtain a non-abelian set of unitary gates. 

The situation is somewhat reminiscent of the "tilted interferometry approach" [25]. There, the 
change of a surface topology with time is used to upgrade topological quantum computation with 
Ising-anyons from non-universal to universal. In our case, the change is from abelian to non-abelian. 
As a final comment, the change of surface topology with 'time' appears as a discontinuous process. 
This is an artifact of the mapping from 3 spatial dimensions to 2 spatial dimensions plus time. In 
the 3D cluster picture there is no discontinuity. 

2.2 Transforming defect configurations 

In the following we discuss equivalence transformations on the defect configuration. Two local 
defect configurations are "equivalent" if they have the same effect in a larger topological circuit. 
The transformation rules allow us to simplify topological circuits and to prove circuit identities. 

The defects are regions in the cluster lattice C but for quantum information processing the 
details of their shape are unimportant. Only the topology of the defect configuration matters. As 
a result, the diagrams of defect strands representing quantum gates such as in Fig. Q] bear a certain 
resemblance to link diagrams. There are indeed similarities but there are differences, too. The 
main similarity is that the line configurations representing defect strands in these diagrams respect 
Reidemeister moves, 



They are valid for both types of defect and all possible combinations. A first difference is implicit 
here: there are two types of lines, primal and dual. 




} dual qubit 



(7) 



} primal qubit 




(8) 
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b) 




1 : ancilla 



2: control In 



3: ancilla + 



4: target In 



control Out 



target Out 



Figure 4: a) Deformed version of the CNOT-gate displayed in Fig. [T^,. b) Equivalent circuit, 
representing a CNOT gate between the control and target qubit. 



Next, we examine the crossings. The crossings of defect strands of the same type are trivial, 



(9) 



Only the crossing of two defects of opposite type is non-trivial; see Eq. ([7]). However, the double 
monodromy of two defect strands of opposite type again is trivial, 



(10) 



There is a special rule for a pair of defect strands supporting a qubit which is encircled by a defect 
of the opposite color. This configuration amounts to measuring the a stabilizer generator or 
([5]), respectively, of the encoded magnetic or electric qubit. This measurement acts as the identity 
operation on the code space , such that 



qubit 



"qubit. 



Q qubit 



quint 



(11) 



So far, it looks as if we were discussing link diagrams with colored components. But there is more 
phenomenology. Three or more defect strands can be joined in a junction. The defect configurations 
thus form graphs. Here is an equivalence transformation by means of which junctions are introduced 
into the configuration, 



^4 




(12) 



This is a somewhat complicated rule. The following happens here: The dual loop on the l.h.s of 
(I12p is contracted. If it has external legs (two are shown), then these are joined in a vertex. The 
primal defect strands passing the dual loop (three are shown) are cut and reconnected. The upper 
and lower parts of each are joined at a vertex. A dual cage is formed around these newly formed 
primal vertices. 

To prove that the two configurations are indeed equivalent it needs to be checked that the set 
of supported correlation surfaces is the same for each. This is beyond the scope of this paper; 
however, one member of this set is displayed in Fig. The equivalence holds for an arbitrary 
number (including none) of involved primal and dual defects. The dual relation (primal defects «-> 
dual defects) also holds. 
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Figure 5: Extended correlation surface passing the junction. The shown surfaces look the same 
far from the location where surgery was performed. If the defect strands pairwise form qubits, the 
shown surface imposes a Z (g) Z-correlation for either defect configuration. 



Finally, simply connected defect regions can be shrunk to a point and removed, 



(13) 



These rules will be used in Section [3] to simplify sub-circuits. To illustrate their use, we give two 
examples of deriving circuit identities in a topological manner. First, A(X) C) t|0) c (0| = It® |0) c (0|. 
In the topological calculus, 
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Second, A(X) ai bA(X)i, ta A(X) a ^ = SWAP(a, b). In the topological calculus, 
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3 Completing the universal set of gates 

The topologically protected gates, the CNOT and preparation/measurement in the X- and Z- 
eigenbasis, are shown in Fig. [TJ The X- and Z-measurements are obtained by reversing the time- 
arrow in the corresponding state preparations. 

We can complete these operations to an universal set by adding exp(i|Z), exp(ijZ), exp(z^X). 
The fault-tolerant realization of these gates requires error- free ancilla states \Y) := (|0) +i\l))/y/2 
and \A) := (|0) + e i7r/4 |l))/\/2. These states are first created in a noisy fashion using the element 
displayed in Fig. [6l and then distilled [30]. For details, see Section [6] and Appendix lAl 
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Once the ancilla states \A) and \Y) have been distilled they are used in the circuits of Fig. [7h,b 
to produce the desired gates. The gate exp^Z) is probabilistic and succeeds with probability 1/2. 
Upon failure, the gate exp(— i^Z) is applied instead, which can be corrected for by a subsequent 
operation exp(ijZ). The latter gate is deterministic modulo Pauli operators, which suffices for the 
QC C . 

Their fault-tolerant QCc-realizations for the above gates are shown in Fig. [7b, d. These realiza- 
tions are obtained from pasting the standard elements for the CNOT and measurement together 
and subsequently applying the defect transformation rules (JEJ) - (|13|) . 



4 Mapping to a two-dimensional system 

The dimensionality of the spatial layout can be reduced by one if the cluster is created slice by slice. 
That is, we convert the 'simulated time'-axis — introduced as a means to explain the connection with 
surface codes — into real time. Under this mapping, cluster qubits located on time-like (space-like) 
edges of C, C become syndrome qubits (code qubits) which are (are not) periodically measured. 

Most important is the region V in which we have topological error protection. Therein, space- 
like oriented A(Z)-gates remain and time-like oriented A(Z) gates are mapped into Hadamard gates. 
The temporal order of operations is displayed in Fig. Note that every qubit is acted upon by 
an operation in every time step. The mapping to the two-dimensional structure has no impact 
on the information processing. In particular, the error correction procedure is still the same as in 
fault-tolerant quantum memory with the toric code. 

In the 3D version, we use |+)-preparations and A(Z)-gates for the creation of \<f>)c-> an d subse- 
quently perform local X, X ± Y, Y and ^-measurements. We now give the complete mapping for 
these operations to the 2+1 dimensional model. 

1. Space-like edges (primal and dual). We group together the respective |+)-preparation, mea- 
surement and trailing time-like oriented A(Z)-gate, and denote the combination by {\+) , A(Z) , P} . 
If the measurement on the trailing end of A(Z) is in the Z-basis, then 

{|+>,A(Z),.P}— P. (14) 

Otherwise, 

{\+),A(Z),P x } — ► H, 
{\+},A(Z),P x±Y } — > He** 2 
{\+),A(Z),P Y } He** 
{\+},A(Z),P z } — P X 



4Z (15) 



2. Time-like edges (primal and dual). For each such edge, we group together the respective 
preparation and measurement, and denote the combination by {|+),P}. Then, 

{\+),Pz} — I, (16) 
{\+),P} — {|+),P}, fovP^P z . UDj 



bare S-qubit 




Out (encoded) 



Figure 6: Preparation of the ancillas \Y) and \A), encoded with the surface code. To obtain \Y) 
and \ A), the singular qubit S is measured in the eigenbasis of Y or (X + Y)/y2, respectively. 
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Figure 7: One-qubit rotations, a) Circuit for performing C/z-gates using the ancillas \Y), \A). b) 
Circuit for performing JTx-gates using the ancilla \Y). c) and d) show the corresponding defect 
configurations. 



3. Space-like oriented A(Z) -gates. 

K{Z) a>h — A(Z) 0j6 . (17) 

Remark 1. No qubit in the scheme is ever idle between preparation and measurement. The 
identity in the first line of Eq. (|16h can be replaced by the 1-qubit completely depolarizing map 
without affecting the scheme. The respective qubit will be re-initialized before its next use. 

Remark 2. From the perspective of information processing, the space-like oriented gates A(Z) a ^ 
in (I17p have no effect if a £ -D V 6 G D. They may consequently be left out. Keeping these redundant 
gates in the scheme, however, does not affect the threshold; see remark 4. We keep the redundant 
A(Z)-gates in order to maintain translational invariance of the (Ising) qubit-qubit interaction. 

Remark 3. For physical realization of the scheme with cold atoms in an optical lattice it may 
be preferable to use a double-layer 2D structure instead of a single layer. The advantage then 
is that all qubits within one layer, including the S-qubits, can be read out simultaneously. One 
'clock cycle' consists of the following steps: 1) Ising interaction/ A(Z)-gates between all pairs of 
nearest neighboring qubits in the lattice; 2) Simultaneous measurement of all qubits in layer a, 
re-preparation of all qubits in layer b; 3) Same as 1); 4) Same as 2), with a <-> b. 

Ideally, one would use a bcc lattice half a cell thick but an sc lattice one cell thick also works. In 
the latter case, some redundant A(Z)gates/ Ising-type interactions and Z-measurements increase 
the number of error sources and thus moderately reduce the error threshold. 



5 Fault-tolerance and threshold 

Error model. We assume the following: 

1. Erroneous operations are modeled by perfect operations preceded/followed by a partially 
depolarizing single- or two-qubit error channel 

T X = (l-p 1 )[I]+p 1 /3([X} + [Y} + [Z}), 

T 2 = (l- P2 )[I]+p 2 /15([X a X b ] + .. + [Z a Z b ]). 

2. The error sources are a) faulty preparation of the individual qubit states |+) (error probabil- 
ity pp), b) erroneous Hadamard-gates (error probability pi), c) erroneous A(Z)-gates (error 
probability p 2 ), and d) imperfect measurement (error probability pm)- 
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a) b) 




Figure 8: a) Temporal order of operations in V after the mapping to 2D. Shown is the elementary 
cell of the 3D lattice with one axis converted to time. The labels on the edges denote the time steps 
at which the corresponding A(Z)-gate is performed. The labels at the syndrome vertices denote 
measurement and (re-)preparation times [£m>*p], and the labels (tn) denote times for Hadamard 
gates. The pattern is periodic in space, and in time with period six. b) The logical cell. It is 
rescaled from the elementary cell of the lattice £ by a factor of A in each direction. The defects 
have cross-sections d x d. 



3. Classical processing is instantaneous. 

When calculating a threshold, we consider all error sources to be equally strong, pi = P2 = Pm = 
Pp := p, such that the noise strength is described by a single parameter p. Storage errors need not 
be considered because no qubit is ever idle between preparation and measurement. 

Error correction. The three relevant facts about fault-tolerance in the QCc are 

• The error correction in V is topological. It can be mapped to the random plaquette 7j2-gauge 
model (RPGM) in three dimensions [21]. If there are non-trivial error cycles of finite smallest 
length I then, below the error threshold, the probability of error etop is 

e top ~ exp(-K(p)l). (18) 

• The topological error correction breaks down near the singular qubits. This results in an 
effective error on the S-qubits that needs to be taken care of by an additional correction 
method. This effective error is local because the S-qubits are far separated from another [5]. 

• The cluster region D need not be present at all. It is initially included to keep the creation 
procedure of the cluster state translation invariant, and subsequently removed by local Z- 
measurement of all qubits in D. The purpose of D is to create non-trivial boundary conditions 
for the remaining cluster. 

The fault-tolerance threshold associated with the RPGM is about 3.2 x 10 -2 [26], for a strictly local 
error model with one source. Also see [27]. The threshold estimates given in this paper are based 
on the minimum weight chain matching algorithm [28] for error correction. This algorithm yields a 
slightly smaller threshold of 2.9% [29] but is compuationally efficient. 

Concerning the exponential decay of error probability in V in Eq. (|18p . the dominant behavior 
is both predicted from a Taylor expansion of etop in terms of the physical error rate p (truncated 
at lowest contributing order [5]) and confirmed by numerical simulation (see Fig. [9]). Beyond the 
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Figure 9: Exponential decay of the failure rate e of the topological error correction as a function of 
the length I of the shortest non-trivial error cycle and the number N of such cycles, at p = 1/3 p c . 



dominant exponential decay there is a polynomial correction, et op ~ exp(— kI)1@. Eq. (36) of [5] 
predicts such a correction and the numerical simulation finds it. However, the exponents (3 differ. 
Eq. (36) of [5] predicts (5 = —1/2 for a strictly local error model. The numerical simulation finds, 
for the close-to-local error model introduced above, [3 = —1.3 ± 0.2 in the time-like direction and 
= —0.9 ± 0.2 in either space-like direction. . Because of the uncertainty in the values of (3 we 
do not include the polynomial correction in our analysis of the operational overhead. This is safe 
because it is to our disadvantage. However, the exponential decay dominates and the effect of the 
polynomial correction is very small. 

The rate k of the dominant exponential decay of error is potentially different along space-like 
and time-like directions, due to the anisotropy of the error model. The numerical simulation finds 
marginal differences at p = p c /3, 



K 
K 



0.85 ±0.03 
0.93 ±0.03 



(time-like) , 
(space-like) . 



(19) 



Error correction in 5*. The <S-qubits are involved in creating noisy ancilla states p ~ |^4)(A|, 
p Y w |y)(y| encoded by the surface code, via the construction displayed in Fig. [6l Due to the 



effective error on the S-qubits, these ancilla states before distillation carry an error eg 
(A\p A \A), e%:=l- (Y\p Y \Y) given by 



1 



6p. 



(20) 



Threshold. There are two types of threshold within the cluster, the topological one in V and 
thresholds from |^4) and |Y)-state distillation in S. An estimate p\ to the topological threshold is 
found in numerical simulation of finite-size lattices, 

p v = 7.5 x 10~ 3 . (21) 

The result of the simulation is displayed in Fig. [TUl 
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Figure 10: Numerical simulation for the topological threshold in V. The curves are best fits taking 
into account finite size effects of the lattice size I. Beyond the smallest lattices, these finite size 
effects quickly vanish, and the curves intersect in a single point to a very good degree of accuracy. 
The value of p at the intersection gives the threshold. 



The recursion relations for state distillation, in the limit of negligible topological error, are to 
lowest contributing order eA_ 1 = 35(e^) 3 (c.f. [30]) and ej +1 = 7(e^) 3 . The corresponding distillation 
thresholds expressed in terms of the physical error rate p are 

P f = ^= « 2.8 x icr 2 , p Y c = ^= « 6.3 x irr 2 . (22) 

The topological threshold is much smaller than the distillation threshold, and therefore the former 
sets the overall threshold for fault-tolerant quantum computation. 

In our previous paper [5], the non-topological threshold was the smaller one because the Reed- 
Muller quantum code was probed in the error correction mode instead of the error detection mode 
associated with state distillation [30]. If we include state-distillation into the setting of [5], the 
fault-tolerance threshold increases to 

p c = 6.7 x 10" 3 , (23) 

which is the topological thresholcd. The threshold (|23|) supersedes the result of [5]. 

Remark 4- The effect of removing the redundant space-like oriented A(Z) gates (c.f. Remark 2) 
is to reduce the effective error on the 5-qubits. Eq. ([20|) then is replaced by = = This 
affects neither the threshold nor the overhead scaling. The distillation threshold increases but it 
already is the larger one. Also, as will be discussed in the next section, the exponent which governs 
the overhead scaling is a geometric quantity unaffected by the values of and in Eq. (I20p . 

6 Overhead 

We are interested in the operational cost per gate, O3, as a function of the circuit size £1. To 
facilitate the calculation of O3 it is helpful to introduce the notions of the scale factor A, the defect 



2 The value ([23} differs from (|21[) due to minor differences in the error model. Specifically, [5] requires a 2D structure 
of three or more layers instead of a single one in the present discussion. Consequently, the used operations and the 
order of operations differ. 
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thickness d, the gate length L, and the gate volume V. 

In the presented scheme, quantum gates are realized by twisting defects. For that purpose alone, 
the defects could be line-like structures. Then, the elementary cell of the lattice C constitutes a 
building block out of which quantum gates and circuits are assembled. However, in such a setting 
the property of error correction is lost due to the presence of short error cycles. To eliminate such 
errors the logical elementary cell is rescaled to a cube of A x A x A elementary cells. The cross-section 
of a defect with the perpendicular plane becomes an area of d x d elementary cells (see Fig. [8b). 

The gate length L is the total length of defect within a gate, measured in units of the length 
of the logical cell. We will subsequently use the gate length for an estimate of the gate error 
remaining after topological error correction. The gate volume V is the number of logical cells that 
a gate occupies, each consisting of A 3 elementary cells of C. Each elementary cell is built with 16 
operations. 

Let £top(G , A, d) be the probability of failure for a gate G, as a function of the scale factor A, 
defect thickness d and of its circuit layout. The operational overhead 03(G) is then 

O a (G) = 16X 3 V G exp (e top (G, A, d)Q) . (24) 

The exponential factor comes from the expected number of repetitions for a circuit composed of Q 
gates G. For a given fi, the overhead should be optimized with respect to choosing X(Q) and d(Q). 

6.1 CSS-gates 

The simplification for CSS gates is that no 5-qubits are involved and all operations are topologically 
protected. To perform the optimization in Eq. (I24p we need to know the gate error etop as a function 
of G, A, d. The errors leading to gate failure may either be cycles wrapping around defects of 
opposite color or relative cycles ending in defects of matching color. The probability of gate failure 
is exponential in the length of the shortest cycle or relative cycle, and proportional to the number 
of such error locations. The minimal cycle length is A(d + 1) and the number of such cycles is 
equal to the gate length XL G . The minimal length of a relative cycle leading to an error is A — d. 
It stretches between two neighboring defect segments one logical cell apart. The number of such 
relative cycles is at most 2L G X(d + 1). There are shorter relative error cycles near junctions, but 
they are homologically equivalent to the identity operation, 

M = b|_ = LL (25) 

Thus, the gate failure rate is 

e t0 p(L G , A, d) = \L G (exp (-4k (d + 1)) + 2(d + 1) exp (-k(A - <£))) . (26) 

We may now use this expression in Eq. (|24j) and optimize for given computational size O. As an 
example, the operational overhead for the CNOT-gate of Fig. [1] is displayed in Fig. [TTJ 



The scaling limit. We now perform the optimization of O3 in (|24p with respect to A, in the limit 
of large circuit sizes fi. First, the gate error etop m (126|) is minimized when both exponentials in f)26f) 
fall off equally fast, i.e. d op t = X op t/5 for large d, A. Further, the overhead O3 in (f24"j) is minimized 
near 

e(A(0)) = I/O. (27) 

Then X op t ~ IuQ/k, and 

ln 3 

Oa ~ — • (28) 
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Figure 11: Operational overhead as a function of the circuit size. The upper curves are for the gates 
exp(f|Z) and exp(z^X), the lower for the CNOT. 



6.2 Non-CSS gates 

The estimation of the overhead for the non-CSS operation is along the same lines but more com- 
plicated, due to the involved magic state distillation. Every level I of distillation is associated with 
its own scale factor A; and a defect thickness d\. The optimization of O3 thus is over the larger set 
of parameters A = { { A , d , Xi , d l , .. . , \ imax ,d ima J,l max }. 

Further, there are now two types of error. First the previously discussed error of topological 
error correction resulting from non-trivial error-cycles far away from the S'-qubits. Second, there is 
error associated with the S-qubits where topological error correction breaks down. 

The distillation of states \A) and \Y) uses S'-qubits at the lowest level. | ^-distillation is based 
on the [15,1,3] Reed-Muller quantum code [30] and the distillation of \Y) on the [7,1,3] Steane 
code. The | A)-distillation is performed using the circuit displayed in Appendix [A] of volume Va 
and length La (See Table [1]). At each level / it requires 15 states \A) of level / — 1 and, on average 
1705/512 3.33 states \Y). It succeeds with a probability of 1 — 15e^ 4 i 1 — et op (LA, A/_i, where 
ef_-^ is the error of the states \A) at level I — 1. If successful, the residual error in the states | A) at level 
I is ef = 35(e^ 4 _ 1 ) 3 + eto P (Ly, A;_i, The distillation for \Y), of volume Vy and length Ly, takes 

7 states \Y) of the next-lower level, and succeeds with a probability of 1 — 7ef_ 1 — et op {Ly, A;_i, 
If it succeeds, the residual error is ej = 7(e^l 1 ) 3 + e top (Ly, A;_i, d[_i). 

The above expressions for success probability and residual errors hold to leading order in the 
contributing error probabilities £-x-\- Further, a gate error cannot simultaneously lead to 

termination of the circuit and to a residual distillation error. Thus, we overestimate both error 
probabilities by adding the full weight e top (L, A/_i , di-i) to them. 

The operational overheads for state distillation at level I, and Oj t , and the corresponding 
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gate 


volume 


length 


CNOT-gate of Fig. QJi (packed) 


V 2 = 12 


L 2 = 22 


Uz-gaie of Fig. [T^ 


V hz = 2 


L\ jZ = 3 


Ux-gate of Fig. [TH 


V 1>x = 4 


Li, x = 4 


|y) -distillation circuit 


V Y = 120 


L Y = 120 


^-distillation circuit of Fig. 1121 


V A = 336 


L A = 362 



Table 1: The gate volume and length for various gates and sub-circuits. 



residual errors ef, ej are described by the recursion relation 



U 3,l 
U 3,l 



1 



l-7e ; y _ 1 -e t0 p(iy,A i _ 1 , 

35(e^ 4 _ 1 ) 3 + e t op(LA, Aj_i,dj_i), 



ro^+mxUvY), 



(29) 



The initial conditions are O A = Ol Q = 16, and (p0|). 

The distillation outputs states |^4) and \Y) at level £ max . One such state |^4) and, on average, 
1/2 state \Y) is used to implement a gate exp(i^Z) via the circuit displayed in Fig. [7^, of volume 
and length -Li j2 . Its overhead is 

OT = (oi Uax + \ol lmax + 24(A,_) 3 V M ) exp (( e f_ + + e top (L M A Z _, d lma J) 0) . 

(30) 

The operational overhead needs to be optimized over the parameter set A. This has been done 
numerically [32] , and the result is shown in Fig. [TT1 



The scaling limit. First, we compare the two contributions to e A in (|29p . 35(e^ 1 ) 3 and et op - 
If etop is much larger than 35(e A _ 1 ) 3 it inhibits the convergence of ancilla distillation. Additional 
distillation rounds are needed which are the most expensive component. If, to the contrary, et op 
becomes much smaller than 35(e^ 4 _ 1 ) 3 it does not help the ancilla distillation anymore but blows 
up the size of the logical cell. Therefore, for optimal operational resources, both contributions 
are comparable. Then, in the large size limit, lne^ 4 = Slne^, A; = 3A/_i, di = 3di-±. Further, 
the success probabilities 1 — 15e^ and 1 — 7eJ for ancilla distillation quickly approach unity with 
increasing distillation level I. Therefore, in the large size limit, for the point of optimal operational 
resources, the recursion relations (|29|) can be replaced by 



0$ 
01 
A 3 
d 
In e A 
\ lne y 



15 





1705 
512 

7 




16 V A 
16 V y 
27 



\ 



V 



3 / 



o A 

A 3 
d 
In e A 
\ lne y 



(31) 



Thus, A p Ol t ~ 27 ? , lnef,lnef ~ 3'. Then, with e ~ I/O Ij27|> . 

O a ,OJ~ (In nf. 



(32) 
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Note that the distillation operations, for the case of perfect CSS-gates, are associated with the more 
favorable scaling exponents log 3 15 ~ 2.46 and log 3 7 ~ 1.77, respectively. However, in our case 
the topological error protection of CSS gates must keep step with the rapidly decreasing error of 
state distillation, by adjusting the scale factor A. This leads to a scaling exponent of 3 for the 
CSS operational resources (c.f. Eq. ([28]) ). which dominates the resource scaling of the entire state 
distillation procedure. 

Discussion. We have found that there is one dominant exponent which governs the scaling of the 
operational overhead for all gates from the universal set, O3 ~ ln 3 r2, c.f. Eqs. (|28p and (|32p . This 
exponent is a geometrical quantity. Its value, 3, derives from the fact that the cluster state used 
in the scheme lives in three spatial dimensions, and that errors are identified with line-like objects 
(1-chains). Details of the implementation such as the volume and length of the distillation circuits 
play no role for the scaling. This summarizes the main results of this section. 

Let us now go beyond scaling and look at the pre- factors. Because of the uniform overhead 
scaling the ratio of operational costs for non-CSS to CSS-gates is constant in the limit of large 
computational size f2. Inspection of Fig. [TT] shows that this ratio is in disfavor of the non-CSS 
gates. 

Without going into much detail, we would like to point out that there is room for improvement 
here. The ratio 0^/0^ X ^ is proportional to Va/cI' 3 . Herein, d! is the shortest length of an error 
that goes undetected in the distillation circuit. Its value is constrained by 1 < d! < d, where d is 
the distance of the used code (3 for the Reed-Muller quantum code and the Steane code. The code 
distance d shall not be confused with the defect thickness d introduced earlier in this section.). In 
the present discussion, c.f. paragraphs preceding Eq. ([29]) . we have used the lower bound d' = 1 to 
simplify the error counting. 

The advise is to replace the Reed-Muller (Steane) quantum code by another [n, k, d] CSS-code for 
which the encoded gates exp(i^Zj) (exp(i^Zj)), i = 1 .. k, are transversal, and with the additional 
properties of having a large distance d and a good ratio k/n. Such codes need to be searched for 
systematically. 

The distillation circuits can be further optimized. The logical circuit depth can be reduced to 
3 for |j4)-distillation and 2 for |y)-distillation, independent of the the code parameters n, k and d. 
The circuit height 3 remains unchanged. Thus, the circuit volumes per output qubit Va, Vy can be 
reduced to V A = 9(n/k + 1), V Y = 6(n/k + 1). 

7 Summary and outlook 

In this paper we have discussed in detail the error threshold and overhead for universal fault-tolerant 
quantum computation based on the one-way quantum computer with a three-dimensional cluster 
state. By conversion of one spatial cluster dimension into time we have reduced the dimensionality of 
the scheme to two. Also, the described scheme only requires translation-invariant nearest-neighbor 
Ising interaction among the qubits. These features should facilitate future implementation. We 
envision cold atoms in optical lattices [8] , two-dimensional ion traps [9] , quantum dot systems [10] 
and arrays of superconducting qubits [TT] as suitable candidate systems for experimental realization. 

On a more abstract level, we have initiated the discussion of the topological properties of the 
defect configurations. We have described a set of transformation rules that allow us to switch 
between equivalent configurations. We have applied these rules to simplify sub-circuits and to 
derive circuit identities, by a sequence of operations reminiscent of the Reidemeister moves for link 
diagrams. 

There are a host of questions that remain open, from the applied to the abstract. Below a few 
are listed. 
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• Optimization of the error threshold. With the current implementation of error correction we 
have exhausted the capabilities of the minimum-weight chain matching algorithm. There is 
one improvement that promises a noticeable gain. So far, error corrections on the mutually 
dual lattices C and C run entirely separate. However, errors on C and C are correlated such 
that error correction could benefit from cross-talk between the two lattices. 

• Transversality of encoded gates. Motivated in part by the discussion of overhead in Section [6] 
but a topic of more general theoretical interest is to find further stabilizer codes that posses 
the capability of performing non-Clifford gates transversally. 

• Robustness of error threshold. Here we have discussed a logic gate-based error model. What 
about more physical error models such as, for example, spins coupled to an Ohmic bath? 

• Connection with the category-theoretic work of Abramsky and Coecke. In this paper we have 
used an encoding with two holes per logical qubit. There is another encoding that gets by 
with a single hole, making additional use of the external system boundary. In that other code, 
for both the primal and the dual defects, the "cups" of Fig. [T|b denote the preparation of a 
Bell state, and the corresponding "caps" Bell measurement. They provide a concrete physical 
realization of the corresponding abstract elements in the category-theoretic calculus of |33j. 
Also, the authors of [33] introduce a "line of information flow" . It is represented by the defect 
strands in our scheme. 

If the teleportation identity was the only phenomenology supported by the defects we would 
not get very far in terms of fault-tolerant quantum computation. To this end, it is crucial to 
have two distinct types of qubits, primal and dual, which interact in a non-trivial manner ([?])■ 
Now, the question is whether this enlarged phenomenology can be included in the category- 
theoretic framework of [33] and whether it enriches that framework. 

Acknowledgments. RR would like to thank Hans-Peter Biichler, Jiannis Pachos, Almut Beige, 
Simon Benjamin, Parsa Bonderson, Trey Porto and David Weiss for discussions. KG is supported 
by DOE Grant No. DE-FG03-92-ER40701. JH is supported by DTO. RR is supported by the 
Government of Canada through NSERC and by the Province of Ontario through MEDT. 

A Circuit for state distillation 

We use a variant of the magic state distillation circuit described in [30], adapted to the QCq- The 
topological circuit is displayed in Fig. [12j The procedure is this: we start with a Bell state, encode 
one of its qubits with the Reed-Muller quantum code and measure each of the 15 qubits leaving the 
encoder in the eigenbasis of Aj — Yj. 

The Reed-Muller quantum code [30] has the nice property that the A-syndrome can be measured 
in the A, Y, A + Y and A - Y-basis. Further, = (g)-^ ^ip- Therefore, through the local 

A — Y-measurements and classical post-processing we can both learn the A-syndrome and project 
the encoded qubit of the Bell pair into an eigenstate of A + Y . Thus, we simultaneously project 
the unencoded qubit into the state X a Z b \A), with a, b £ {0,1} depending on the measurement 
outcomes and on which of the four Bell states was used. We keep this qubit if the above A- 
syndrome measurements yield a trivial outcome. In this case, the residual error q is, to leading 
order, ef = 35(e^ 1 ) 3 (c.f. [30]). 

The local A — Y-measurements are performed by a unitary operation exp(— i^Zi) followed 
by an Admeasurement. Each such unitary requires one ancilla \A) and, with probability 1/2, 
one additional ancilla |Y) such that one round of magic state distillation performed in this way 
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consumes 15 states \A) and, on average, 15/2 states \Y). With a small modification^!, we can reduce 
the average number of required |Y)-states to 1705/512. 

The distillation circuit for |Y)-states is constructed in a similar manner. It is based on the 
Steane code and requires seven input states \Y) in each round. 
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Figure 12: QCc-realization of the |^4)-state distillation (a). The dots on the defect lines are the 
ports to connect \A)- and \Y) states for Z-rotations (c.f. Fig. [7]). The main part in the distillation 
is the encoder for the [15, 1, 3] Reed-Muller quantum code, displayed as a quantum circuit in (b). 



B Effective error model on C and C 

After the mapping described in Section [J] the physical setting is in two dimensions. However, the 
topological error correction is still performed on the three-dimensional lattices C and C. The error 
model of Section [5j including gate error for one and two-qubit gates, preparation and measurement, 
effectively results in Z-errors on individual edges and correlated errors on two edges of C or C, 
respectively. The location of correlated errors is shown in Fig. [T3l 
The effective error channel on time-like edges is 

r 1>t = {{i-u 2 )[i t ] + y 2 [z t )) o2 o{{i-i PP )[i t ] + i PP [z t ))o 

((i-y M )[i t \ + l PM [z t ]). 

The effective error channel on a space-like edge — horizontal or vertical — is 

?i.< = ((l-^2)[/ s ] + ^ 2 [^])° 3 o((l-|pi)[/ s ] + |pi[^])° 2 . (34) 

3 Denote by J the set of subsets of {1,2, ..,15} such that ® ieJ Xi is an encoded gate (including the identity 
operation) for all J £ J. Then, for the Reed-Muller quantum code, the Clifford unitary 7 exp(i-| Zi) is also 
an encoded gate VJ £ J, namely I or exp(— i\Z\ Now suppose that after probabilistic implementation of the 
7r/8-phase gates, 7r/4-phase gates on a set K C {1, 2, .., 15} are required. Then, it is equivalent to perform Uk®j — 
<S> i€K ex P(*f z i) <8>i 6 j exp(if Zi), VJ G J, modulo local Pauli operators Zi. (Note that exp(if Z)\A) = X\A) = \A).) 
We minimize the support of Uk®j mod {Zi} by varying J 6 J. In this way, we reduce the average number of 
|Y)-states required in a distillation step to 1705/512 « 3.33. 
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horizontal vertical time-like 

Figure 13: The correlated errors (same for on the primal and dual lattice). Shown are horizontal, 
vertical and time- like faces of C and C. Thick lines indicate error locations. 



The effective error channel for each of the correlated errors displayed in Fig. [13] is 

r 2 = + 



(35) 



The probability of time- like and space- like individual errors is different. In the error correction 
procedure this is accounted for by using non-uniform weights in the minimum- weight chain matching 
algorithm [28]. Likewise, the correlated errors in the boundary of space-like faces are accounted for 
by including additional diagonal edges. 
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